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Abstract 

European options can be priced when returns follow a Student's t-distribution, 
provided that the asset is capped in value or the distribution is truncated. 
We call pricing of options using a log Student's i-distribution a Gosset ap- 
proach, in honour of W.S. Gosset. In this paper, we compare the greeks 
for Gosset and Black-Scholes formulae and we discuss implementation. The 
i-distribution requires a shape parameter v to match the "fat tails" of the 
observed returns. For large v, the Gosset and Black-Scholes formulae are 
equivalent. The Gosset formulae removes the requirement that the volatility 
be known, and in this sense can be viewed as an extension of the Black- 
Scholes formula. 
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1 Introduction 



It is known that Student's t-distribut ions fit the returns of stocks and eg- 



uity indices better than normal distributions ( 


Mandelbrotlll963 


. iFama 


1965. 


Praetz 




1972. Blattbere and Gonedes 


1974. 


de Jong and Huismanll2000 


. Platen and Sidorowicz 


20071. IZhu and Galbraith 


2001 


}. Gerier et al. 120091. IWu and Kerchevall 12010). 


Prices 



for European options can be calculated for returns that follow a Student's t 
distribution i f the value of the underlying asset is capped or if the t-distribution 



is truncated ( Bouchaud and Sornettd Il994j . iMcCauley et al.\ 120071 . ICassidy et al. 



20091 ). Capping the value of the asset or truncating the distribution keeps the 



integrals, which are required to price the option, finite. 

Truncation is not an academic exercise. It is unphy sical for the value 
of an asset to approach infinity fjBouchaud and Sornettd Il994l ). Thus truncation 
is required for a model of values of assets and of returns on assets to reflect re- 
ality. Capping has a similar effect as truncation and is simply a slightly different 
approach to creating a physical model for the price of an option. 

It is the "fat tails" of the Student's ^-distribution that make trunca- 



tion necessary ( iBouchaud and Sornettdll994l . lMcCaulev et q/.ll2007l ICassidv et al. 



20091 ). The normal pdf, on which the Black-Scholes formula for the price of an 



European option is based, decreases rapidly for values that are more than sev- 
eral standard deviations beyond the mean value. This rapid decrease keeps the 
integrals, which are required to price an option, finite and thus allows a price to 
be found without consideration of truncation or capping. However, this rapid de- 
crease in the tails is not consistent with the observed distribution of returns, and 
thus calculations based on a normal pdf underestimate the probability of events 
that are several standard deviations beyond the mean. 

The Student's ^-distribution offers support from — oo to +oo (lEvans et al. 



1993). However there is no requirement that the full domain be used when working 
with a distribution. If the data do not or can not exist over the full domain and 
yet the frequency of occurrence of the data do follow the distribution over a subset 
of the full domain, it is correct to truncate the domain to match the observations 
or physical constraints. One need only use a concept from conditional probability 
to maintain the normalization of the truncated pdf. In terms of the stock market, 
the magnitude of the returns and the value of an asset are truncated or capped. 
In practice, one does not observe a stock that is worth an infinite amount or daily 
returns on indices that are > 30%. 

In working with a Student's ^-distribution to price European options, 
it is necessary to decide upon the shape parameter v and the scale parameter 
<jt (the scale parameter ot is the volatility) for the t-distribution, and to decide 
where to place the cap or the truncation. The decision on the value for the cap 
or truncation is based to a large extent on tolerance to risk or equivalently on 
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the anticipated maximum value for the asset. One can choose to work with, e.g., 
p = 99.999% confidence or 95% confidence, or choose the corresponding maximum 
value x c {p). The shape parameter v and the volatility a t are a different story. 
These parameters should be chosen to match the underlying distribution. 

In this paper we discuss pricing options when the underlying distribution 
is a Student's i-distribution. We call a model to price European options when 
the underlying distribution is a t-distribution a Gosset model, in honour of W.S. 
Gosset who published under the name Student f Zabelll 20081 ). We briefly discuss 
implementation, present the Gosset equations in Appendix B, discuss what is an 
appropriate shape parameter v based on historical returns, and compare some of 
the greeks for the Gosset and Black-Scholes prices for European call options. 



2 Implementation 

The Gosset formulae for the prices of European options are similar to the 
Black-Scholes formula. In the limit as the shape parameter v approaches infin- 
ity, the Gosset formulae become the Black-Scholes formula. Both formulae are 
based on tabulated functions. The Black-Scholes formula requires values for the 
error function. The Gosset formulae require tabulated values for integrals of log 
Student's ^-distributions. Fortunately both the error function and the values for 
integrals can be calculated numerically in short time periods. 

Appendix A is a copy of the output of a Maple worksheet that was 
used to determine the time required to calculate and print to the screen the value 
of a call option for 100 different values for the asset S. This script took 20 s to 
execute. One hundred evaluations of the Maple intrinsic function blackscholes(S', 
K, r, T, a) for 100 different values of S would take 5 ms to complete and to 
print to the screen. The execution time for pricing a call with a Gosset formula 
could be reduced if optimized code is used to evaluate the integrals, similar to the 
optimized code that is used for evaluation of the error function to price an option 
using the Black-Scholes formula. 

One challenge with use of the Gosset formula is that spreadsheets do 
not have built-in functions to evaluate the integrals, unlike the functions that can 
be used to price an option with the Black-Scholes formula. However, in time this 
shortcoming can be remedied. 

One advantage that the Gosset formulae enjoy over the Black-Scholes 
formula is the inclusion of the "fat tails" that returns display. The Gosset formulae 
can be viewed as extensions to the Black-Scholes formula. The Gosset formulae 
price European options where the returns are normally distributed but the true 
(parent) volatility is not known a priori (i.e., the true volatility is estimated as 
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the sample standard deviation) or the true volatility is a random variable and 
is distributed as an inverse chi variate. The Black-Scholes formula is based on 
a volatility that is known; there is no uncertainty in the volatility in the Black- 
Scholes approach. It is known t hat it is important to include sto c hastic volatility 
in pricing options - see , e.g., (IHestonl Il993l . iBaskhi et ail 119971 . iMoriconil 120071 . 



Wu and Kerche vail 1201 Of ). 



For the Gosset formulae to hold strictly, the volatility must be dis- 
tributed as a chi random variable or as an inverse chi variate. The requirement of a 
chi or inverse chi pdf depends on the approach that is cho sen to justify t heoretically 



the f act that returns follow a Student's i-distribution (IPraetzl I1972L IGerig et al 



20091 ). One can postulate that the Student's t-distribution and the concomitant 



fat tails arise from the uncertainty introduced by estimating the true volatility 
from the sample standard deviation, or one can postulate that the volatility is 
a random variable. Presumably the assumption that the volatility is a random 
variable takes into account the uncertainty in estimating the true volatility from 
the sample standard deviation. 

Historical data sh ow that the yolatility follows, to a good approximation, 
an inverse chi distribution (IPraetzl Il972l IGerig et a/.ll2009l ). which suggests that 
the true volatility is a random variable. The volatility is not constant in time. 
The markets show periods of calm and of wild swings in prices. These periods of 
calm and wild swings in price can not be predicted with certainty and thus the 
volatility is a random variable to an observer of the market. 

We find that both the inverse chi and the chi distributions fit the volatil- 
ity data for the S&P 500 and the DJIA, with the inverse chi distribution providing 
a moderately better description of the volatility than the chi distribution. Thus 
one would expect the returns to be described by a Student's t-distribution and not 
a normal pdf, as is observed. Hence the predictions of the Gosset formulae should 
be preferred over the predictions of the Black-Scholes formulae for the prices of 
European options. 

Figure 1 gives the price of call and put options as a function of the 
strike price K, asset value So = 50.00, confidence level p = 0.999, risk free rate 
r, time to expiration T, rxT = 0.03, volatility Of = 0.3, and various values 
of the shape parameter v for the Student's t-distribution. The prices for v = 
oo were calculated with the Black-Scholes formula. The solid curves for v = 
3, 8, and 21 were calculated using the Gosset formulae, Eqs. (B10) and (Bll) of 
Appendix B, for a truncated pdf. The broken lines were calculated using the Gosset 
formulae, Eqs. (B6) and (B7) of Appendix B, for a capped distribution. Clearly 
the difference between the Gosset price and the Black-Scholes price increases as 
v becomes smaller. As v becomes smaller, the tails of a Student's t-distribution 
become "fatter". It is also clear from Fig. 1 that the capped Gosset price, i.e., 
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Figure 1: Plot of the option price as a function of the strike K. 



the price for an option that is calculated using a value for an asset that is capped 
and follows a log Student's t-distribution, is marginally greater than the truncated 
Gosset price. The additional cost for a capped value allows for the possibility that 
the value of the asset may exceed the cap value. Truncation assumes that the 
value of the asset will not exceed the cap. 

Table 1 gives the critical values x c (p) and increase in the value of the 
asset at the critical value for p = 0.999 for four different values of the shape 
parameter v. The critical value x c and confidence level p are defined such that the 
probability of drawing a return < x c is p. 

To price an option with either the Black-Scholes approach or the Gosset ap- 
proach, the value of the asset is assumed to vary as exp(<7i£), i.e., the logarithm 
of the value of the asset is distributed as £. Table 1 shows that truncating the 
Student's t-distribution at p = 0.999 with cry = 0.3 allows for a maximum increase 
of the value of the asset of 21.421 times for v = 3. For the normal distribution v = 
oo and at x c the relative increase in the value of the asset is 2.527 times. Table 1 
also shows that the truncation need not be restrictive, and thus pricing of options 
with the Gosset formula is realistic. It is worth noting that one could truncate or 
cap at different values of p, say p = 0.99999, which would allow for a maximum 
increase of the value of the asset of 106 times for v = 3, a? = 0.3, or say p = 
0.95, which would allow for a maximum increase of the value of the asset of 2.075 
times for z/ = 3,crr = 0.3 The choice of the level for the cap or the truncation 
might depend on the underlying asset. Penny stocks might increase dramatically 
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Table 1: Critical values and maximum increase in the value of an asset for p = 
0.999 and cry = 0.3 for various shape parameters. 





v = 3 


v = 8 


v = 21 


V = oo 


critical value, x c 


10.215 


4.501 


3.527 


3.090 


exp(0.3xx c ) 


21.421 


3.858 


2.881 


2.527 



in terms of percent increase, but little in terms of absolute value, and thus require 
a large critical value to capture the possible changes. 



3 Estimation of the Shape Parameter v 



One question that arises naturally in pricing an option with a Gosset for- 
mula is what value to use for the shape parameter v. Under the assumption that 
the returns are normally distributed, v < N where N is the number of indepen- 
dent samples that have been used to estimate the volatility. The argument for the 
maximum value is straightforward. 

Let s be the sample standard deviation that has been calculated for N 
independent samples from a normally distributed parent population. The sample 
standard deviation s follows a chi distribution with v = N — 1 degrees of freedom , 
Xu(s), for data that are normally distributed (jEvans et a/.lll993l . lKreysziglll972l ). 



Any measurement of s has an uncertainty and this uncertainty in s is reflected 
in the shape parameter of the sampling distribution. For normally distributed 
samples where the standard deviation of the parent distribution is estimated from 
the sample population, the sampling distribution is a Student's ^-distribution with 
v degrees of freedom, f s (x, v). 

The uncertainty in the estimation of the true (parent) standard devia- 
tion can be significant. 

The (sample) volatility is the sample standard deviation for the daily 
retur ns, and a symm etric 95% confidence interval for the true (parent) volatility 
a is riKrevszidll972h : 0.77xs < a < 1.43xs for v = 21; 0.83xs < a < 1.27xs for 
v = 43; and, 0.87xs < a < 1.17xs for v = 87, where s is the sample standard 
deviation (volatility) with v degrees of freedom that has been calculated from the 
daily returns and it has been assumed that the returns are normally distributed. 

The maximum value for u, assuming that the volatility is estimated 
as the sample standard devia tion of the re turns, is obtained when the volatility 
is constant in time. Praetz (jPraetzl Il972[ ) and Gerig et al. (IGerig et all 120091 ) 
found that the variance of returns is a random variable that changes slowly in 
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time. Our simple analysis o f historical re turns, which we pr esent below, agree s 



with the findings of Praetz (IPraetzl Il972l ) and Gerig et al. ( iGerig et al\ 120091 ) . 



This randomness in the volatility adds weight to the tails of the distribution of the 
returns. 

Additional uncertainty in the volatility leads to a lower value (i.e., fatter 
tails for the underlying distribution) for the shape parameter in the Gosset formula. 
The additional uncertainty arises when one wishes to price an option for exercise at 
some time T in the future. Since the volatility is a random variable, the volatility 
at time T, a?, is unknown and the additional uncertainty must be taken into 
account. 

We estimate values for the shape parameter v from historical data from 
the DJIA (Oct 1928 to Feb 2009) and S&P 500 (Jan 1950 to Feb 2009) equity 
indexes. 

In Tables 2 to 5 are mean and median values for 22 day volatilities. 
For the simulations, 22 samples were generated and the volatility (i.e., the sample 
standard deviation of the 22 samples) was calculated. The process was repeated 
1000 times, and the mean, median value, and the standard deviation s of the 1000 
22 day volatilities were calculated and are reported in the tables. 

For each group of 22 values, the minimum value and the maximum value 
were dropped to give 20 values, and the volatility was calculated for this truncated 
set of data. The process was repeated 1000 times and the mean, median value, 
and the standard deviation of the 1000 volatilities were calculated, as for the full 
set of 22 values. 

For each group of 20 values, the minimum value and maximum value 
were dropped to give 18 values. The calculations described above were repeated 
and the numbers entered into the tables. 

For calculation of the expected volatility, the pdf was assumed to be 
truncated at a critical value x c such that the probability of an observation x in the 
range — x c < x < x c was equal to N/22, i.e., P{— x c < x < x c } = pn = N/22. 
There likely is a rigorous method to calculate the standard deviation for dropping 
the minimum and the maximum values from a finite number of random draws from 
a parent population, but the method that was adopted was simple and shows the 
trend. 

The expected volatility (a) was calculated as 



{a) = r< with Pn = r< m d ^ (1) 

V J-xc Pn J-xc 

where /(£) is the appropriate probability density function and the denominator 
Pn maintains the normalization of the truncated pdf. 

Table 2 gives results of the simulation and calculation for the procedures 
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Table 2: Simulated medians and standard deviations of 22 day volatilities for 
samples drawn from a t-distribution with v = 3 and from a normal pdf with 
variance = 3, and the expected volatility (a). 





simulation 


calculation (a) 


N 


t, v — 3 


normal, a 2 =3 


t, v = 3 


normal, <r 2 =3 




median 


s 


median 


s 




<*> 


22 


1.48 


0.65 


1.70 


0.28 


1.73 


1.73 


20 (drop 2) 


1.07 


0.25 


1.41 


0.25 


1.01 


1.39 


18 (drop 4) 


0.87 


0.19 


1.21 


0.24 


0.816 


1.18 



Table 3: Simulated mean, standard deviation s, and median of 22 day volatilities 
for samples drawn from a ^-distribution with v — 3 and from a normal pdf with 
variance = 3. 





£ (simulated), v — 3 


normal (simulated), a 2 =3 


N 


mean 


s 


median 


mean 


s 


median 


22 


1.61 


0.65 


1.48 


1.71 


0.27 


1.70 


20 (drop 2) 


1.09 


0.25 


1.07 


1.43 


0.24 


1.41 


18 (drop 4) 


0.88 


0.19 


0.87 


1.22 


0.23 


1.21 



described above. Note that the normal distribution has been scaled to have the 
same variance as the ^-distribution. Note also that the median 22 day volatility de- 
creases as the maximum and minimum samples are dropped, and that the median 
value for the ^-distribution decreases faster than the median value for the normal 
distribution as the maximum and minimum samples are dropped. The standard 
deviation of the 1000 22-day volatilities decreases greatly for the ^-distribution 
as the maximum and minimum samples are dropped, but changes little for the 
normal distribution. This owes to the fat tails of the ^-distribution with v — 3. 

Table 3 shows the means, the medians, and the standard deviations for 
1000 samples for the same ^-distribution and normal distribution in Table 2. The 
means are in agreement with the expected values in Table 2 and show the same 
progression as the median when the minimum and maximum values are removed 
from each sample of 22 data points. 

The data in the tables show that the fat tails of the ^-distribution result 
in large uncertainties of the volatility (i.e., the sample standard deviation s of the 
22 day volatility is large for the ^-distribution). The fat tails are also evident in 
the difference between the median and the mean values for N = 22. Note that 
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Table 4: Comparison of the normalized S&P 500 mean and median 22 day volatil- 
ities (676 values) to the normalized expected volatilities. 





S&P 500 


expected volatility (a) 


N 


mean 


median 


t, v = 3 


normal 


22 


1 


1 


1 


1 


20 (drop 2) 


0.803 ± 0.010 


0.81 


0.584 


0.803 


18 (drop 4) 


0.677 ± 0.011 


0.68 


0.471 


0.683 



Table 5: Comparison of the normalized DJIA 22 day volatilities (918 values) to 
simulations. 





DJIA 


simulations 












normal (mean 


N 


mean 


median 


t, V 


= 3 


or median) 


22 


1 


1 


mean 


median 


1 


20 (drop 2) 


0.799 ± 0.008 


0.802 


0.67 


0.72 


0.83 


18 (drop 4) 


0.675 ± 0.009 


0.680 


0.55 


| 0.59 


0.71 



the uncertainty in the volatility, as measured by the standard deviation s, reduced 
dramatically for the data from a ^-distribution after truncation of the minimum 
and maximum values, but barely changed for a normal distribution. 

Table 4 compares values of the volatility for the S&P 500 daily returns 
to the simulations. The S&P data was broken into 676 non-overlapping segments of 
22 values. The values were sorted and then the volatility (i.e., the sample standard 
deviation) was calculated for the middle 22, 20, and 18 values of the 22 returns. 
The mean, median, and standard deviation s of the volatilities were calculated 
over the 676 volatilities. For ease of comparison, each column was normalized by 
the value for N = 22. 

Table 5 compares means and medians of the 22 day volatilities for the 
DJIA equity index to simulations. 



Tables 4 and 5 show that the short-term DJIA and S&P 500 volatility 
behaves more like a normal pdf than a Student's ^-distribution. The interpretation 
is that for a short term (22 returns), the returns are approximately normally 
distributed. However, the volatility changes with time, and for longer periods of 
time, the uncertainty in the volatility comes into play and the distribution becomes 
less normal like and more t like. This is evident from the data in Table 6, where 
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Table 6: Normalized mean and median volatilities for the DJIA (459 and 230 
values) and the S&P 500 (338 and 169 values). 





DJIA 


S&P 


N 


mean 


median 


mean 


median 


44 


1 


1 


1 


1 


40 


0.773 ± 0.012 


0.792 


0.776 ± 0.016 


0.793 


36 


0.646 ± 0.012 


0.658 


0.649 ± 0.016 


0.648 


88 


1 


1 


1 


1 


80 


0.748 ± 0.017 


0.748 


0.753 ± 0.023 


0.786 


72 


0.622 ± 0.016 


0.625 


0.626 ± 0.021 


0.663 



the normalized means and medians are reported for calculations over longer time 
frames of 44 and 88 returns. 

Table 6 shows that the means and medians of the truncated samples are 
smaller for longer time periods. This suggests that the data are departing from 
normally distributed. 

An error propag ation equation was used to est imate the uncertainty 
in the ratio of the means ( Bevington and Robinsor] 2003 ). The uncertainty in 



B/A was calculated as 2xsb/a, where (sb/a) 2 = (sb) 2 /A 2 + (sa) 2 (B / A 2 ) 2 — 
(B / A 3 )sab 2 , to give a 95% confidence level assuming that A and B are normally 
distributed. The covariance s^ B was in the range of 0.97 to 0.91 for A equal to the 
mean for the full range and B equal to the mean for dropping the bottom and top 
5% or 10% of the samples used to form A. 

Figure 2 is a plot of the normalized expected volatility as a function of 
the shape parameter v for the Student's t-distribution. The long tic marks at 0.803 
and 0.683 are the normalized expected volatilities for samples that are normally 
distributed. The upper curve is for = 20/22 = 0.9090 while the lower curve is 
for pjv = 18/22 = 0.8182 The curves are normalized by the expected volatility for 
p N = 22/22 = 1.000 

The normalized mean volatilities for the S&P 500 and DJIA equity 
indexes were placed on Fig. 2 to match the normalized mean value of the volatility 
with the normalized expected volatility curves. This placement on the expected 
curve provides a means to estimate the shape parameter and the uncertainty in 
the shape parameter. The normalized mean volatilities are found in Tables 4, 
5, and 6. The vertical bars give the 95% confidence intervals for the normalized 
mean values. The horizontal bars give the uncertainty in the shape parameters 
for the uncertainty in the mean value. From Fig. 2 it appears that a Student's 
^-distribution with a shape parameter v — 13 ± 4 would be appropriate for a 44 
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day volatility and that a Student's t-distribution with a shape parameter u = 8 ± 
2 would be appropriate for an 88 day volatility, based on the historical data used 
in the analysis. For a 22 day volatility, the shape parameter should be < 21. 

The data on Fig. 2 at v = 2.65 ± 0.11 and v = 3.24 ± 0.19 were 
obtained from fits to the full data sets for the DJIA and S&P 500 equity indexes 
(jCassidy et a/.ll2009l ). The shape parameters were obtained by fitting frequency 
of occurrence data to a Student's t-distribution. For Fig. 2 the returns were 
sorted from smallest to largest and the volatility sa was calculated as the standard 
deviation over the N = 20186 returns in the DJIA data set and N = 14870 returns 
in the S&P 500 data set. The bottom 4.545% and top 4.545% of the data were 
truncated, and the standard deviations sb were calculated. Normalized standard 
deviations of 0.583 ± 0.006 and 0.634 ± 0.007 were obtained by forming sbAa 
and are plotted on the figure. The bottom and top 4.545% returns were truncated 
from the previously truncated lists, and the normalized standard deviations were 




Figure 2: Standard deviation of the truncated data as a function of the number of 
degrees of freedom v. 
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obtained. The uncertainties in the normalized standard deviations were estimated 
as 95% symmetric confid ence intervals f or a chi-squared variate with degrees of 
freedom equal to N — 1 (jKrevszigl If 9721 ). The normalized data obtained in this 
manner, where the degrees of freedom were obtained from independent fits to 
frequency of occurrence data, are consistent with the expected volatility curve. 

Given the above, it appears that the foll owing model d escribes well the 
returns, in agreem ent with the results of Praetz fjPraetzl Il972l ) and Gerigei al. 
fjGerig et a/.ll2009l ). The returns are normally distributed with a defined volatility 
for a short term (of order 22 consecutive returns). The volatility is a random 
variable that changes slowly in time. A i-distribution, which is a mixture of a 
normal pdf and a chi distribution for the reciprocal of the standard deviation (i.e., 
1/volatility) describes the pdf for the returns now and at some time in the future, 
where the volatility is uncertain owing to sampling and to the slow and random 
development of the volatility in time. In mathematical notation, the sampling 
distribution f s (x, v) = ff(x\a) Xvi^ 1 ) da is a Student's t-distribution when the 
parent distribution f(x\a) is a normal distribution with standard deviation a and 
cr follows an inver se chi distribution (or, equivale ntly, the reciprocal of a follows a 
chi distribution) flPraetzlll972l . iGerig et a/il2009h . 

Figure 3 is a plot of a fit of a chi distribution to the reciprocal of the 
22 day volatilities for the S&P 500 equity index from Jan 1950 to Feb 2009. It 
can be observed that the reciprocal of the historical 22 day volatilities are, for the 
most part, distributed as a chi variate. Since the returns over 22 days are found 
to be normal like (see the tables) with a given volatility but the reciprocal of this 
volatility is distributed as a chi varia te, the distrib ution of the returns is expected 
to follow a Student's t- distribution (|Praetzlll972r) . A Student's ^-distribution for 
the returns is observed (jCassidy et a/.ll2009i ). In Fig. 3 the data deviate slightly 
from the best fit curve in the tail on the right hand side. The upward sloping 
curves are the cumulative density functions (CDF) for the data (red) and the best 
fit function (black). For a perfect fit, the black and red CDF curves would overlap. 

Figure 4 is a plot of a fit of a chi distribution to the 22 day volatilities 
(the volatilities, not the reciprocals of the volatilities) for the same S&P 500 data 
shown in Fig. 3. It is interesting to note that the chi distribution fits the 22 day 
volatilities reasonably well. The fit deviates in the tails on the right hand side, 
more so than the fit of Fig 3. Fits of a chi distribution to synthetic data from an 
inverse chi distribution show that the inverse chi data has a fatter tail on the right 
hand side, similar to the fat tail on the right hand side that the S&P 500 data 
show in the figure. 

The shape parameter for the i-distribution is u, the number of degrees 
of freedom. From the analysis presented here, the shape parameter is related to 
estimations of the distribution of the volatility. If the volatility does not change 
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0.0 100 200 300 400 

1/(22 day volatility) 

Figure 3: Data and fit of a chi distribution to the frequency of occurrence of the 
reciprocal of the 22 day volatilities for the S&P 500. 




S&P 500 



0.0 0.01 0.02 0.03 

22 day volatility 

Figure 4: Data and fit of a chi distribution to the frequency of occurrence of the 
22 day volatilities for the S&P 500. 
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in time and is known (i.e., the volatility is not estimated as the sample standard 
deviations of the returns), then v = oo and one recovers the Black-Scholes formula 
for the price of an European option. The Black-Scholes formula is derived on the 
assumption that the volatility is constant and known. If the true (parent) standard 
deviation (volatility) is not known and the returns are normally distributed, then 
a Student's ^-distribution should be used. If v is small, then it is understood that 
the true volatility might take one of a broad range of values. 

The simple analysis presented also shows what might be reasonable 
values (see Fig. 2) for the shape parameter of the Student's t-distribution. Since 
the volatility is estimated as the sample standard deviation of the daily returns, 
it is expected that v < N , where N is the size of the sample that was used to 
compute the volatility. 



4 Greeks for Truncated Student's ^-Distributions 
and Capped Values 

The Gosset formulae for the price of an option by capping the value of 
the asset or by truncating the pdf are given in Appendix B. These formulae are 
extended to include a floor or a truncation at the floor level. The notation is 
defined in Appendix B. 

The expression for Delta for a truncated Student's ^-distribution, t A, is similar 
to the expression for A for the Black-Scholes formula. In the limit of no truncation 
(Pp = 0, p c = 1, x c = oo) and v — oo, the expression below reduces to A for the 
Black-Scholes formula. 

t Q Sq jMtZ,K T ^S )-rxT ^ x _ K ) 

For the Black-Scholes formula, t Z = exp(|o"T 2 ), Pc — Pp = 1, £ c = oo, and the 
integral can be written in terms of the error function. 
The expression for C A for a capped asset price is 



A_ 9cCo _ r efrrQ/ r (0# e^'°>x(l-p c ) 

- q JlnUZxKr /S )-r X T C Z C Z U 



C 



i) .s. 



Figure 5 plots t A for Pp = 0,p c = 0.999, K T = $49.00, rxT = 0.03, and 
ot = 0.3 as a function of So and for v — 3, 5, and 21. For comparison, A for the 
Black-Scholes formula is shown. On the scale of the figure, t A for v > 20 lies on 
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top of the Black-Scholes A. The dashed lines plot Delta for a capped distribution, 
C A. 

t A is a function of p c . For large v, t A shows a weak dependence for 
values of p c near unity. This is expected since for large v there is little difference 
between the Black-Scholes and Gosset formula. The Black-Scholes formula is in- 
sensitive to values in the tails and thus one could reasonably expect the Gosset 
formula for large v to also be insensitive to values in the tails. Figure 6 plots t A 
for v = 3, p p = 0, K T = $49.00, rxT = 0.03, and ot = 0.3 as a function of S 
and for p c = 0.99, 0.999, and 0.9999. A small value of v was chosen to highlight 
the differences. The dashed lines plot Delta for a capped distribution. 

tA is also a function of the volatility ot- Figure 7 plots t A for u = 3, 
p p = 0, K T = $49.00, rxT = 0.03, and p c = 0.999 as a function of So for o T = 
0.1, 0.2, and 0.3 The dashed lines give C A. 

An expression for t r is given below. This expression reduces to the 
Black-Scholes formula for no truncation and v = oo. 



d \C K T xe-™ T xf r ( 



ln( t ZxK T /S )-rxT 



(4) 

" dS A S z xa T x(p c -p p ) KJ 

Figure 8 plots t F for the same parameters as used to create Fig. 5. The 
value of t r decreases as the shape parameter decreases. 4 r is also a function of 
p c and a T . Figure 9 plots t T for v = 3, p p = 0, p c = 0.999, K T = $49.00, rx T 
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Figure 6: Delta as a function of the value of the asset S for different levels of 
truncation p. 
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Figure 7: Delta as a function of the value of the asset S for different volatilities. 
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Figure 8: Gamma as a function of the value of the asset S for different values of 
the number of the degrees of freedom v. 



= 0.03 as a function of Sq for o~t = 0.1, 0.2, and 0.3 The dashed lines gives c r, 
which is T for a capped asset. The expression for c r is 

_ d\C KT X e-* r X f ( McZ»K T M-r»T } 

c - dS 2 - S 2 xa T (5) 

which is similar to the expression for t T. 

An expression for Vega for a truncated t-distribution, t V, is given below. 

tV = 7^7 = ,Zx(p!-pJ Kl t ZK T /S )-rxT Ce^^/r-CO d£ 

ln( t Z K T /S )~rxT C ^MOdt (6) 

Figure 10 is composed of plots of Veg Qb clS 8b function of Sq for the standard 
parameters used in this work. The expression for Vega has a term owing to the Z 
in the limits of integration. This term can not be neglected. Vega increases as the 
shape parameter for the t-distribution decreases. The curve for v = 40 lies just 
above the curve for the Black-Scholes Vega. The dashed lines are for Vega for a 
capped distribution. An equation for Vega for a capped asset, C V , is 
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9 ffj- 



i( c ZA' T /S )-rxT 

J ln( C ZK T /S )-rxT 

, g oXe ^T^) x (i_ Pc) T / _ i 9eZ 



-Z 2 9 or 



(7) 



ftf = Pp x x p e^^) + jj / r (£) ^ + (1 _ Pc ) x x c e("**«) 



The derivatives are messy. The number of terms and parameters makes 
it difficult to determine by looking at the expressions the importance of terms and 
parameters. As a practical matter, it is probably easier to compute the derivatives 
numerically as [ t C (Ci> C2, ••• Ci+e, Ci+i, •••) - tC ((i, C2, ••• Cu d+i, ■■■) ]fe fc> r 
a range of parameters and plot the values to build an understanding. 

Figure 11 is plots of G as calculated numerically as a function of S'o for 
the standard parameters. Theta increases as the shape parameter of the Student's 
^-distribution decreases. The dashed lines are for a capped distribution. Theta 
was calculated as C (p, S, K, rx (366/365), a T x (366/365) - 5 , v) - C (p, S, K, 
r, cr T , u). 

Figure 12 is plots of the numerical derivative of the price of a call with 
respect to the shape parameter v as a function of S'o for the standard parameters. 
This derivative is zero for an option priced with the Black-Scholes formula. Not 
surprisingly, the price of a call is sensitive to the shape parameter for small v. The 
dashed lines are for a capped distribution. 
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Figure 11: Theta as a function of the value of the asset S for different values of 
the number of the degrees of freedom v. 
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Figure 12: Derivative of the price of a call option with respect to the number of 
degrees of freedom v as a function of the value of the asset S for different values 
of the number of the degrees of freedom. 



Figure 13 is plots of the logarithm of the derivative of the price of a 
call with respect to the truncation p = p c as a function of So for the standard 
parameters. The dashed lines are for a capped distribution. Note the large value 
of the magnitude of the derivative for p approaching unity. Note also that the 
maximum value for e is << (1 — p) and that the change in the price of a call for 
a change in p of e is the derivative xe. 

The integrals that are required to price a call option contain terms 
with exp(<7r x c ). As p c approaches unity, x c approaches infinity, and the integrals 
required to price a call option assuming the returns follow a Student's t-distribution 
approach infinity. The case of x c approaching infinity is unphysical as this suggests 
that the value of the asset exp(ax x oo) must also approach infinity. As p c 
approaches unity, the price of the option increases to account for the values of the 
asset that are being included. 



5 Conclusion 

It is known that a Student's t-distribution fits the returns of stocks better 
than a normal pdf. A Student's i-distribution has "fat tails", and the shape 
parameter v and scale parameter <jt of the t-distribution can be selected to match 



20 



4 - 



CD 

£ 3 

CO 

> 

CD 
T3 



CD 
O 



1 1 1 1 1 1 1 1 1 








v = 3,p= 0.9999 


— v = 3, p = 


0.999 








v = 3, p = 0.99 — 


v = 21,p 


= 0.999 




1 1 


1 1 1 


1 1 1 1 



2 - 



1 - 



42 



46 



50 
S 



54 



58 



Figure 13: Plots of the logarithm of the derivative of the value of a call option 
with respect to the level of truncation p as a function of the value of the asset S 
for different values for the level of truncation. 



the fat tails of the returns. 

Prices for European options can be determined for assets that follow a 
log Student's t-distribution if the value of the asset is capped or if the ^-distribution 
is truncated. Capping or truncating the distribution keeps the integrals that are 
required to price an option finite. We call the formulae that use a Student's t- 
distribution to calculate the values of options Gosset formulae, in honour of W.S. 
Gosset. Gosset used the pseudonym Student to publish his work. 

In the limit of large v the Gosset formulae converge to the Black-Scholes 
formula for the values of European options. The Gosset formula can be viewed 
as an extension of the Black-Scholes formula to allow for a volatility cx^ that has 
uncertainty or is a random variable. The Black-Scholes formula is based on the 
assumption that the volatility is known and constant. The returns must be nor- 
mally distributed and the volatility must follow a chi or inverse chi distribution for 
the Student's t-distribution to hold strictly. The historical returns from the DJIA 
and S&P 500 appear to follow these distributions as the Student's t-distribution 
fits well the returns on these indexes. 

It is necessary to know the shape parameter v to use the Gosset formu- 
lae. 

The maximum value of v < N where N is the number of independent 
returns that are used to calculate the volatility. The shape parameter, in general, 
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will be less than the maximum value for large N (say TV > 22). The volatility is 
not constant in time, and for large N, the uncertainty of the volatility broadens 
the pdf that describes the returns and hence reduces v. Also, for longer time to 
expiration T, the volatility at time T is known with less accuracy, and a reduced 
v is necessary to capture the uncertainty in the future value of the volatility. Fits 
to the returns from the DJIA (Oct. 1928 to Feb. 2009) and the S&P 500 (Jan. 
1950 to Feb. 2009) give v of order 3. Thus it appears that typically v should 
be in the range of ~3 to A 7 "— 1. A method to estimate an appropriate value of 
v by comparison of the standard deviations of the returns as the minimum and 
maximum values are dropped is given. 

Some of the greeks for the Gosset and Black-Scholes formulae for pric- 
ing European options are given and compared. The differences between the two 
approaches are the largest for small v. This is not surprising, given that the dif- 
ference between the two underlying distributions for the returns are the greatest 
for small v. 

New greeks are required to characterize the price sensitivity to param- 
eter change for the Gosset approach. The Gosset formula requires a cap or trun- 
cation and a shape parameter v. Changes in the value of a call option for changes 
in v and the truncation or cap level are presented. 
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A Maple script to price a call option for 100 different values of the asset price is 
presented in this Appendix. The script was used to find the execution time that 
is reported in Sec. 2, Implementation. The intrinsic function to price a call with 
the Black-Scholes formula is included as a comment in the loop. 



> sigma := 0.30: # volatility 
r := 0.03: # risk free rate 
K := 49.00: # strike 

nu := 21: # shape parameter 

p := 0.9999: # confidence level 

_prepare to price with log Student's t 

> eq := proc(nu,xi) 

GAMMA ( (nu+1) /2) /GAMMA (nu/2) /sqrt(Pi*nu) / (l+xi*2/nu) A ( (nu+1) /2) 
end proc: # define procedure to evaluate Student's t-distribution 
eq('nu',t); # display the function 



> t_call := proc (p,xc,S,K,r, sigma, nu) local Z, a, b; # use CAPPED pdf 
b := exp (sigma*xc) * (1-p) ; # evalf( Int() ) = numerical integration 
Z := evalf (Int (exp (sigma*xi) *eq (nu,xi) ,xi = -inf inity . . xc) ) + b; # 
a := (In (Z*K/S) -r) /sigma; # lower limit of integration 

( S/Z*evalf (Int ( exp (sigma*xi) *eq (nu,xi) ,xi = a..xc)) 

- exp(-r) *K*int(eq(nu,xi) , xi = a..xc) + S*b/Z - exp (-r) *K* (1-p) ) ; # 
end proc: # end of procedure to price a call 

> xc := f solve (int (eq(nu, t) , t=0 . . x) +0 . 5 - p, x) ; # find critical value 
for S from 1.00 to 100.00 do # evaluate 100 times 

C := t_call (p,xc, S ,K, r, sigma, nu) ; # evaluate using Student's t 
# C := blackscholes (S , K, r, 1.0, sigma): # intrinsic function 
printf ("%6.2f\t%6.3f\n" ,S,C) # print result 



> restart: # comments in green and follow a # symbol 




end do : 



xc :=4.492860131 



Figure 14: Maple script for Appendix A. 
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B 



The equations to price an option by capping the value of an asset or 
truncating the underlying pdf are presented in a general form. The approach 
follows the work of Cassidy, Hamp, and Ouyed ( jCassidy et a/.ll2009l ). 



B.l Pricing the Option 

Let St be the price of a stock at time t, t > 0. Let be the strike price 
at time T, where T is the time when the option expires. 

Let St = A t exjp((Tt£) be the value of the stock where £ is a random 
variable and let the probability density function for the return be / r (£). If £ is a 
Student's t variate, then 

/r(0de= r((V;i)/g x _!?L 
M?; r>/2)Vi^ + e 

The value of a European call option, calculated at the time of expiration T, 
is Ct = E{(5V — Kt) + }i which is the expectation of the maximum value of {St 
— Kt, 0}. The value of a European put option is Pt = E{(Kt — 5V) + }. These 
expre ssions for the values of the options follow from the arbitrage theorem ( Ross! 



20071 ). The values of the options at time t = are obtained from the expected time 
value of money. If r(t) is the risk free rate, then Cq = E{CV x exp(— Jr(t)dt)} 
= Ct x exp(— rx T) when the risk free rate is assumed to be time independent, 
with a similar equation for Pt- 

The average value of a stock at time t = when one wishes to price 
an option is Sq. The average value of the stock at some time later is EjS^} 
= E{A t exp(a t £)} where the drift is contained in A t . For a martingale (a fair 
wager), EjS^} = S , xexp(rx t). One may use the Doob decomposition to offset 
the drift and require E{A t exp(a t £)} — E{A' t exp(a t (£— £ ))} = Soxexp(rxt). This 
is equivalent to requiring that the pdf for St should be centred about £ D , the drift 
owing to the risk premium. £ Q is not a random variable. The exp(rxi) term takes 
into account the time value of money and allows for comparison of values at two 
different points in time. 

The difficulty with pricing an option for which the return follows a Stu- 
dent's t-distribution is that E{5t} is infinite. The exp(a t £) term in the expectation 
dominates over the behaviour of the pdf for large £. Two methods for pricing 
an option are presented below. One method is to place a cap on the value of the 
asset. The other method is to truncate the pdf. Truncation is physical because 
it is not possible for the price of an asset to approach infinity. The equations in 
this appendix are written in general terms. The approach to price an option by 
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capping or truncating is not restricted to a Student's ^-distribution. The approach 
can be used with other probability density functions. In this appendix, provision 
is made for a floor x p or for a truncation at the floor. 

Let p p = P{£ < x p } and p c = P{£ < x c } where x c and x p are critical 
values. In terms of the pdf, p p and p c are equal to (with i = p or c) 

Pi= H fAOdt (9) 

J — oo 

For Pp = 0, p c = 1, and / r (£) a standard normal variate, the equations 
for the prices of European options presented reduce to the Black-Scholes formula. 
We call the formulae for prices of options that have been obtaine d with frit) & 



Student's t-distribution Gosset formulae, in honour of W. S. Gosset (IZabellll2008f ). 



B.2 Cap and Floor 

Assume a floor St = x p for St < x p and a cap St = x c for St > x c and let 
C Z = p p x e {uTX v ] + f X " e {aTi x / r (f) x d£ + (1 - p c ) x e (<TTXc) (10) 

J Xp 

A pre-subscript of c or t indicates that a quantity has been calculated with 
a cap or a truncation. One could write Z as f c Z to indicate a floor as well as a 
cap. However, the extra symbol is not necessary here and is omitted. 
For a fair wager, the expected value of St is found as 

E{S T } = S xe^= c A T x c Z (11) 
and the unknown c At can be found 

C v JrxT) 

C A T = SoXe 7 (12) 
The price of a European call option for a capped asset value is 
cC T = m KT/c A T) ( c A T e^V _ Kt ) x /r(e) x rfe 
+ (l-pl) x ( ^e^-ifr) 

assuming that ln(i^ r / c y4 T )/o"r > x p . 

The price of a European put for a capped asset value is 



(13) 



.Br = P P x (K T - C A T e ^ x ^ ) 

ln(K T / C A T ) Q4] 

+ U ° T (K T - c A T e^^) x / r (0 x d£ 
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In the limit as p p tends to zero, the floor disappears. Note that a cap 
or truncation is necessary to price an option with a Student's ^-distribution. The 
cap or truncation keeps the integrals, which are needed to price the options, finite. 
The floor is not necessary but may make for realistic pricing. The price of an 
option is quantized by the currency. For values less than the smallest unit of the 
currency, it may be appropriate to set the price equal to zero, i.e., to use a floor. 

Put-call parity holds for truncated or capped options. The equation for 
the price of a put is presented for completeness. 



B.3 Truncation 

Assume that S T only takes values x p < S T < x c and thus / r (£) = if £ < 
x p or £ > x c . Let 

tZ = J- ^M|M (15) 

Jx p p c ~ Pp 

The denominator follows from the definition of conditional probability and is re- 
quired to maintain the normalization of the pdf after the truncation: P{x < £ < x+dx \ x p < £ < x c } 
= P{x < £ < x+dx fl x p < £ < x c }/P{x p < £ < x c } = friO^x/iPc — P P ) if Xp < 
£ < x c and zero otherwise. 

The average value of the stock at T is E{S T } = S' xexp(rxt) = t A T 

x t Z or 

t A T = (16) 

The price of a European call option for a truncated pdf is 
fXc (tA T e^-K T )xf r ^)xd^ 

tCT ~ v~^p (17) 

assuming that \n(K T / t A T )/a T > x p . 

The price of a European put for a truncated pdf is 



ln( X tAT) ( K t -tA T e^^)x / r (£) x d£ 



t P T = / ^ V ^ (18) 

Jx p p c — Pp 

Note that t Ax is known in terms of So, r, and tZ(p p (x p ), p c (x c )) and thus 
the option prices can be determined. Note also that t Z ^ C Z . 
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